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Abstract 

We examine how the Sunyaev-Zel'dovich (SZ) effect of clusters can probe Cl . Using 
self-consistent models of X-ray clusters in the context of hierarchical models of structure 
formation, we show that both the mean Compton y parameter and the number of clusters 
observed via the SZ effect strongly depend on Cl . These quantities are higher in open 
cosmological models due to the earlier epoch of structure formation. Specifically, we com- 
pare two models which are able to reproduce the present abundance of X-ray clusters, one 
for = 0.2 and the other = 1. For = 0.2, < y > exceeds the current FIRAS limit 
of 1.5 X 10~ 5 . Also, the SZ redshift source distribution considerably differs in the two 
cosmological models. Thus, these results show that both an improvement in our knowl- 
edge of CMB spectral distortions and the realization of millimeter surveys to look for SZ 
sources over a large area of the sky should produce interesting constraints on and on 
the evolution of the baryonic fraction of virialized structures. 
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1 Introduction 



In open cosmological models, structure formation starts earlier, with its rate approaching a 
constant for redshifts zj ~ l/0 — 1. Thus, the structure observed today already existed at 
Zf. An interesting feature is that structure formation rates are independent of the specified 
spectrum [12]. So, cluster distribution functions should present this caracteristic fi dependence. 
Clusters of galaxies are known to have hot intra-cluster gas which emits in X-rays and is the 
main reservoir of the baryons that fell into the cluster potential wells. This gas will produce 
a unique spectral distortion in the CMB spectrum by its interaction with CMB photons by 
inverse Compton scattering [18]. 

The change, for a given frequency and for a given line of sight, in the sky brightness relative 
to the CMB intensity, which we will call the SZ surface brightness of a source, is given by 
S v = yj l/ (x) } where j v (x) describes the frequency dependence and the Compton y parameter 
measures the distortion magnitude. It has been shown that the integrated SZ surface brightness 
of a source (the total increment or decrement in the normal CMB flux) depends linearly on the 
total hot gas content of the cluster and therefore only on their baryonic fraction [2] and not on 
the spatial distribution of the cluster hot gas : 

S v = (8^jh^)fAx)fBn^Mti 3 [^]^(l + z)D-'(z), (1) 

where z, and Mi 5 are the baryonic fraction, redshift and mass of the cluster. D a and f v (x) 
are respectively the dimensionless parts of the angular distance and of the characteristic SZ 
frequency function. We have assumed h=0.5 in all subsequent calculations. 



2 X-ray Data and the Compton y parameter 

In hierarchical models like the cold dark matter (CDM) family of models of structure formation, 
it is natural to adopt the Press and Schechter (PS) mass function [17] : 

-JlnM = J-±-v(M,z)( - 71 ^ e~" /2 dlnM, (2) 
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where p is the the present mass density of the Universe, v = 6 c D g (z) / 'er(M) , where 6 c D g (z) and 
cr(M) are, respectively, the linear density contrast for collapse at an epoch z and the amplitude 
of the mass fluctuations. The mass variance <r(M) is expressed in the Fourier plane by : 

1 POO 

<r(M) = — / P(k)WUkR)k 2 dk. (3) 
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Here, Wf is some window function and M is the mass contained in a sphere of radius R. In the 
case of a CDM-like perturbation spectrum, the original primordial power-law spectrum (P p ) is 
modified by the presence of cold, non-baryonic dark matter, among other factors. The change 
in the spectrum is contained in the CDM transfer function, for which P(k) = T 2 P p (k) } where 
P p oc k just describes the scalar invariant primordial perturbations. A good fitting formula to 
T(k) is provided by equation(4), where its shape is controled only by the so called T factor [3] 

T(q) = %( ( 1 2 + 4 2 g) 34g) [l + 3.89, + (16.1,) 2 + (5.46,) 3 + (6.71,) 4 ] , (4) 
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Figure 1: The temperature function for the two models; in both cases we have assumed h = 0.5. For each £1 
the curves correspond to the best \ 2 fits to b and T. The local power spectrum index is given by n. 



where q = k/(hT). In the CDM family, T is usually identified as fi/j, but here is considered as 
a free parameter. Also, for the spectrum normalization, we take the bias b = l/a(M) at 8/i -1 
Mpc . 

Because clusters are virialized structures, the mass-temperature relation derived in the frame- 
work of the spherical collapse model, T oc Af 2 / 3 (1 + z) [7], together with the mass function 
yields the temperature function dn/dT. To compare the X-ray data [12], [9], to the theoretical 
temperature functions for different cosmological models, a % 2 fit was performed assuming as 
free parameters the normalization b and the shape of the power spectrum T. 
Figure 1 shows the temperature function dn/dT, with the best % 2 fits to the data for the fi = 1 
and fi = 0.2 cosmological models analysed here. An analogue method was used by Oukbir et 
al. [13], [14] who made a similar % 2 fit to X-ray data, considering the perturbation spectrum 
around cluster scales as a power law. They chose as free parameters b and the local power spec- 
tral index n. The values found by both analyses agree very well, with the b — Y fit indicating 
the same local spectrum index found by Oukbir et al. 

The mean global Compton distortion, < y >, will be the sum of all single distortions pro- 
duced by each cluster. Naturaly, it will be defined as an integral over the entire mass distribution 
function (in this case the PS mass function) : 

< y >= f B / dz— ca T {l + z) 3 / dlnM— — ± —tr 1 , (5) 
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where x 1S the number of electrons per baryon. For X-ray observations indicate values 
as high as 0.05h~ 3 / 2 [6], [8], [19], which agree with nucleosynthesis for fi = 0.2 — 0.3, but 
are inconsistent for fi = f [20]. Also, note that gas in low mass structures is able to cool in 
less than one Hubble time, so the mass integral in equation(5) has a low mass cut-off. This 
cut-off is fO 12 Af and is almost independent of redshift [5]. Typically, at redhifts higher than 
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Figure 2: b — T plane. Thin lines represent y isocontours. The gray zone is the forbidden region : < y > is 
higher than the current FIRAS limit of 1.5 x 10 -5 . Thick contours represent the la and the 2a contours of the 
X 2 fit to temperature function data if errors were gaussian. Asterisks are the best fit values. Note the change 
in the bias scale between the two panels. 



z ~ 5, Compton cooling is efficient, so gas within the structures will cool and, consequently, 
the contribution at these high redshifts to < y > will be small. 

Figure 2 shows X-ray and < y > constraints in the b — T plane. The ellipses (thicker contours) 
represent the lcr and 2<r contours of the % 2 fits. The bias is quite well constrained, unlike T. 
This is not surprising, because b is defined at the cluster scale, while T is a shape parameter 
linking all scales. The best fit values for T are quite different from the value of ~ 0.25 founded 
by the galaxy distribution [f6] but still largely compatible within the fcr confidence level. It 
is clear that for fi = 0.2, only a small region of the b — T region allowed by X-ray data does 
not violate the FIRAS limit of 1.5 x 10~ 5 [15]. Further improvement in the limit on < y > (a 
factor of 2) could provide stronger constraints on this model. 



3 Source Counts 

Integrating the PS mass function, we can study the redshift distribution function of sources 
brighter on the sky than a certain threshold, dN/dft/dz(> SV), and the the corresponding SZ 
source counts. The source counts are expressed by : 

dN , ri . f , dV f llir dn 

— > S v ) = / dz—— / d\l — . 6 

ail J dzdil JM mm (s V) z) dM 

where M m i n is given by equation(f). Figure 3 presents the calculated redshift distribution and 
SZ source counts at 0.75 mm. The panel on the left shows how different the redshift distribution 
of the two cosmological models can be : in an open universe, the long high-redshift tail reveals 
the geometrical effects in the rate of structure formation absent in the flat model. This effect 
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Figure 3: The SZ source counts and the redshift distribution for several S v . The parameters b and T used are 
the best fit values to dn/dT. The calculation is for a wavelength of 0.75 mm. 



is independent of cluster evolution and it is also visble in X-ray redshift temperature function 
analysis, [11], [13]. The righthand panel shows the SZ source counts. A difference of 10 is visible 
at 50 mJy, the expected threshold for COBRAS/SAMBA. If we use a more conservative T value 
of 0.25, the difference in the two distributions is still higher than a factor of 4, thus showing 
the fi„ effects. 



4 Conclusion 

Unlike X-rays, the SZ effect is insensitive to the radial distribution of hot gas in clusters; it de- 
pends only on the cluster baryonic fraction. Thus, the Sunyaev-Zel'dovich effect is a good probe 
of the Universe's geometry. Both the < y > and the SZ source counts show large differences 
between the flat and open cosmological models considered here. It is clear that an = 0.2 uni- 
verse which reproduces the temperature cluster function is left with a small region in the b — T 
plane. This suggests that an improvement in the limit on < y > should be very constraining for 
open models. The redshht distribution of SZ sources is an elegant method of probing fi , due to 
its independence on the normalisation of the perturbation spectrum. Also, as hrst pointed out 
by Korolev et al. [10], cluster source counts should also be a good probe of fi , even if they are 
more sensitive to evolution than the SZ redshift distribution dN/dftdz. In fact, the construc- 
tion of a catalogue of SZ sources could reveal the fi geometrical effects as presented in Figure 3. 
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